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The strain induced pseudo-magnetic field in supported graphene deposited on top of a nanos- 
tructured substrate is investigated by using atomistic simulations. Step, elongated trench, one 
dimensional barrier, spherical bubbles, Gaussian bump and Gaussian depression are considered as 
support structures for graphene. From the obtained optimum configurations we found very strong 
induced pseudo-magnetic fields which can reach up to ~ 1000 T due to the strain-induced deforma- 
tions in the supported graphene. Different magnetic confinements with controllable geometries are 
found by tuning the pattern of the substrate. The resulting induced magnetic fields for graphene on 
top of a step, barrier and trench are calculated. In contrast to the step and trench the middle part 
of graphene on top of a barrier has zero pseudo-magnetic field. This study provides a theoretical 
background for designing magnetic structures in graphene by nanostructuring substrates. We found 
that altering the radial symmetry of the deformation, changes the six-fold symmetry of the induced 
pseudo-magnetic field. 



I. INTRODUCTION 

In most of the experiments on graphene, the 2D atomic 
layer is placed on top of a substrate, which at atomic scale 
is not flat. Geometrically structured substrates affects 
various properties of graphene p], Q, and can prevent 
the crumpling of graphene which occurs for free standing 
graphene without support [3] . Recently, the modification 
of the properties of graphene on top of a substrate were 
investigated. It was found that substrates can induce 
corrugations .modify the electric conductance and deform 
graphene 0,111 . 

Tomori et al used pillars made of a dielectric mate- 
rial placed on top of a substrate which is then over- 
layed with graphene to generate non-uniform strain on 
a micro-scale [6|. Elastic deformations in graphene cre- 
ates a pseudo-magnetic field which acts on graphene's 
massless charge carriers The resulting variation 

of the hopping energies can be viewed as an induced 
pseudo-magnetic field which enters in the Dirac equation. 
Engineering the right topology of the induced pseudo- 
magnetic field can provide magnetic confinement which 
confines electrons in specific regions in space [l(| [H| • It 
has been shown theoretically that inhomogeneous mag- 
netic fields are able to confine massless Dirac fermions in 
a monolayer graphene sheet [12]. 

Here, we investigate several nano-structured substrates 
with different geometrical deformations. We carried out 
molecular dynamics simulations at T=300 K to minimize 
the energy and find the optimum profile of the deposited 
graphene on top of different nanostructured substrates. 
Elongated trench, barrier, bubble, Gaussian bump and 
Gaussian depression are considered as examples of nano- 
structured substrates. The adhesion of the substrate to 
the deposited graphene can induce a very strong pseudo- 
magnetic field which we found depends on the imposed 
boundary conditions on the graphene sheet. Strong 
pseudo-magnetic fields (~ 1000 T) are found around the 



deformed regions in graphene. A substrate with: (i) a 
step, forms two magnetic-barriers around the step with 
opposite sign, (ii) a trench forms two narrow magnetic- 
barriers around the trench boundaries with the same sign 
and one with opposite sign within the trench, and (iii) 
a one dimensional barrier forms two pairs of magnetic- 
barriers around the barrier's wall. The magnetic confine- 
ment for a Gaussian depression in the substrate looses the 
six-fold symmetry of the pseudo-magnetic field which is 
not the case for GE on top of a Gaussian bump. 

This paper is organized as follows. In Sec. II the details 
of the atomistic model are presented. In Sec. Ill we 
present the strain induced gauge field model. In Sec. IV 
we present results for the gauge fields and the pseudo- 
magnetic fields, for various nano-structured substrates. 
The results are summarized in Sec. V. 



II. ATOMISTIC MODEL 

In order to find the optimum configuration of graphene 
(GE) on top of various nanostructured substrates we em- 
ployed classical atomistic molecular dynamics simulation 
(MD). The second generation of Brenner's bond-order 
potential [13| is employed for carbon-carbon interac- 
tion and the van der Waals (vdW) interaction between 
GE and different substrates is modeled by employing the 
Lennard- Jones (LJ) potential, i.e. 



u(r) = 4e[(a/r) 12 - (a/r) 



(1) 



where r is the distance between the two particles, e 
and a are the 'energy parameter' and the 'length pa- 
rameter 1 , respectively. To model the interaction be- 
tween two different types of atoms such as the carbon 
atom (C) and the substrate atom (S), we adjust the 
LJ parameters using the equations ct = ^J~e~cl and 
(Jt = (crc + cr)/2. For carbon we use the parameters 
<jc =3.369 A and ec =2.63 meV. For the substrate atoms 



we set a =3.5 A and e =10.0 meV, which is typical e.g. 
for a Si02 substrate [l4|. The simulation is done for a GE 
sheet with dimension l x =19.17 nm and l y =19.67 nm at 
T=300K. The number of substrate atoms is M =6000. 
In order to model the substrate, a (100) surface having 
a typical lattice parameter £=3 A is assumed. The den- 
sity of sites in the substrate is = £~ 2 - The details 
of the found deformations are reported in our previous 
study [H. 
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FIG. 1: (Color online) The optimum configuration of arm- 
chair graphene over a step located at x = with supported 
longitudinal ends. The colors indicate the size of strain. 



III. STRAIN INDUCED PSEUDO-MAGNETIC 
FIELD 

Generalizing the Dirac equation, which governs the low 
energy electronics of graphene, to curved surfaces is an 
interesting development which may model some cosmo- 
logical problems [8, 9]. The metric of the curved surface 
enters now into the Dirac equation. The origin of the de- 
formations are external stresses which deform graphene 
so that the nearest neighbor distances become non-equal. 
Notice that the external stresses can be induced by the 
substrate. The latter results in modified hopping param- 
eters introduced in the tight-binding model which are 
now a function of the atomic positions t(r) [l6j . Assum- 
ing small atomic displacements (i.e. u = r • — < ao 
where ao is the carbon-carbon bond length) and rewrit- 
ing the Dirac Hamiltonian in the effective mass approx- 
imation with nonequal hopping parameters tells us that 
the strain induces an effective gauge field 

A — — (Uxx ^°yyi C ^H j xy)' l (2) 

where j3 (~2-3) is a constant and u a p is the strain tensor 
including out of plane displacements [8] . The correspond- 
ing pseudo- magnetic field perpendicular to the x—y plane 
is obtained as 

B = d y A x -d x A y . (3) 

This is the pseudo-magnetic field which the electron ex- 
periences in the K valley. We will find B by mak- 
ing the necessary differentiations numerically for longi- 
tudinally supported boundary conditions. Here we are 
mostly interested in the out-of-plane contributions of the 
pseudo-magnetic field which mainly appears around the 
deformed parts of graphene. The other in-plane terms 
contribute less to the pseudo-magnetic field around the 
deformed parts, particularly when the system is larger 
than the size of these deformed parts and is supported 
from boundaries. Notice that in order to perform the nu- 
merical differentiations (in Eq. (2) and Eq. (3)) one needs 
a reference graphene lattice (r^) in order to compare the 
optimized lattice (r-) with the reference system. We used 
the optimized graphene profile at the given temperature 
over a flat substrate as the reference system. However, 
when the boundaries are free there is considerable differ- 
ence (at the boundaries and for some particular systems) 



between the optimized graphene over the deformed sub- 
strate and the reference system. This is due to the fact 
that at finite temperature the free edges of graphene over 
the substrate can vibrate and deform ( due to the substrate 
induced strain) freely while they will not be deformed 
in the reference system. Therefore the reference system 
with free boundaries for some of the systems can be very 
different from the optimized graphene over the deformed 
substrate at finite temperature, hence the differentiation 
is not well defined. Therefore, in this paper we focused 
will be on systems with fixed boundaries, which were 
studied in our previous paper [15], where we have a true 
reference system suitable for numerical differentiations. 

IV. RESULTS AND DISCUSSION 

In this study we investigate several different geometries 
for the substrate which can be realized experimentally. 
For all studied cases we first obtained the optimum con- 
figuration of GE on top of the different nanostructured 
substrates using MD simulations (those results were pre- 
sented in our previous work [15]). Then, for the sup- 
ported boundary condition, we calculate the correspond- 
ing gauge field from which we obtain the pseudo-magnetic 
field. 



A. Step 

An interesting substrate configuration is a step which 
was recently studied in an experiment to measure the 
electronic and morphology of deposited graphene [17j 

h s (x,y) = h o 0(x), (4) 

where Q(x) is the Heaviside step function and ho =1 nm is 
the height of the step. GE with arm-chair direction is put 
on top of the step. In Fig. [I] the optimum configuration 
of GE along the arm-chair direction with longitudinally 
supported boundary condition is shown when placed over 
a sharp step defined by Eq. (j4]). 

The induced gauge field as obtained from Eq. (j2j) is av- 
eraged over the ^/-direction and is shown in Fig. [2fa). All 
atoms at the step region are stretched which results in 
considerable gauge fields around x « 0. Fig. [2fb) shows 
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TABLE I: A list of all relevant parameters used in the paper. 



e,cr 
\,6(x) 
R 
h 

hi , d 
iz a/3 ,A,B 



The graphene length and width 

The energy and length parameters in the van der Waals (vdW) potential for the substrate atoms, Eq. (pQ) 

The wave length and the step function 

The radius of the Gaussian bump or depression 

The amplitude of sinusoidal waves or height (depth) of Gaussian bump/bubble/barrier (depression or trench) 
A shift or vertical distance between graphene and substrate and the width of the trench/barrier 
Strain tensor, strain induced gauge field and magnetic field 
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FIG. 2: (Color online) The averaged gauge fields (a) and 
the induced pseudo-magnetic fields (b) averaged over the y- 
direction for graphene on top of a step as shown in Fig. [T] 
which has been supported from the longitudinal ends while it 
can freely moves along the z-direction. 



the averaged pseudo-magnetic field over the ^-direction, 
versus x. In order to calculate averages we made a 
histogram where l x is divided into 60 equal parts. Notice 
that the induced pseudo-magnetic field is mostly concen- 
trated beyond x=0 and consists of a positive and an ad- 
jacent negative barrier with total average zero. Because 
of thermal fluctuations (i.e. T=300K) the positive and 
negative barrier are only approximately identical. The 
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FIG. 3: (Color online) The optimum configuration of arm- 
chair graphene over a trench located at |x| < 1.5 nm where 
both longitudinal ends were supported in the x — y plane. The 
colors indicate the size of the strain. 



larger the curvature the larger the magnetic field. The 
large pseudo-magnetic field around the step separates the 
GE sheet into a left and a right hand side, where 'B' is 
small. Electrons will be trapped in this region into snake 
orbits and electrons passing perpendicular to this rectan- 
gular part will experience large pseudo-magnetic fields. 
Notice that by changing the height of the step (ho), we 
are able to control the size of the magnetic barrier and 
consequently the magnetic confinement. 



B. Trench 

The other important substrate that we study here is 
an elongated trench 

h s (x,y) = h 9(x 2 -d 2 ), (5) 

with two walls of height 1 nm located at x = ±d = 
±1.5 nm. In Fig. [3] we show the optimum configuration of 
arm-chair graphene with supported boundary condition 
on top of the trench defined by Eq. ([5]). 

The absolute value of the induced gauge field as ob- 
tained from Eq. (j2j) is averaged over the ^-direction and is 
shown in Fig. IHa). All atoms at both sides are stretched 
toward the well region which results in a considerable 
gauge field around x ~ ±d. Fig. HJb) shows the y- 
averaged pseudo-magnetic field. Notice that there is a 
non-zero (|A|) and (B) within the trench which is a con- 
sequence of the bent (non-flat) graphene in the middle re- 
gion. Notice that the pseudo-magnetic fields are smaller 
than those obtained for a step profile. Indeed supporting 
GE longitudinally from two ends prevents GE to move 
into the well and consequently there will be less varia- 
tions in the heights. The magnetic filed profile consists 
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FIG. 4: (Color online) The averaged gauge field (a) and the in- 
duced pseudo-magnetic field (b) averaged over the ^/-direction 
for graphene on top of a well as shown in Fig. [3] which has 
been supported from the longitudinal ends while it can freely 
move along the z-direction. 



of a positive B-barrier inside the trench and two negative 
barriers located at the steps. The total average magnetic 
filed is also zero in this case. 



C. Barriers 

A barrier in the middle of the substrate is reverse sit- 
uation of the previous case. An elongated barrier in the 
^/-direction is parameterized as 

h s (x,y) = h o 0(x 2 -d 2 ), (6) 

with two walls at x = ±d = ±1.5 nm of height of 1 nm. 



FIG. 5: (Color online) The optimum configuration of arm- 
chair graphene over an elongated cubic barrier with \x\ < 
1.5 nm where the zig-zag edges were supported in the x — y 
plane. The colors indicate the size of strain. 



Figure [5] shows the optimum configuration of arm-chair 
GE in the case of supported boundary condition over the 
barrier. 

The induced gauge field as obtained from Eq. ((2]) was 
averaged over the ^-direction and is shown in Fig. [6fa). 
All atoms at both sides are stretched towards the bar- 
rier region which causes considerable gauge fields around 
x « ±d. Fig. EKb) shows the averaged pseudo-magnetic 
field over the ^/-direction which is less than 20 Tesla. Both 
gauge and pseudo-magnetic fields are comparable with 
those found for the substrate with a single step placed in 
the middle of the GE sheet. The main difference is the 
formation of a zero-magnetic field channel in the region 
\x\ < d. The electrons will be trapped in this rectan- 
gular channel which can be realized in experiments. On 
both sides of this magnetic channel there are two double 
magnetic barriers of similar shape. Because of thermal 
fluctuations the barriers are not identical. 



D. Spherical bubble 

The next important deformation of the substrate that 
has been realized experimentally [13, LL9| is a bubble 



(7) 



where R is the radius of the bubble and p 2 = x 2 ± y 2 is 
the radial distance of the i th atom from the center. In 
order to create an uniform discrete atomistic structure 
for the bubble, we set hi = —R/2 where R =2nm. The 
optimum configuration for the longitudinally supported 
graphene over the bubble substrate is shown in Fig. [71 
Due to the supported end GE is elongated longitudinally 
along the supported direction, see inset in Fig. [3 

In Fig. [HJa) the induced gauge field (corresponding to 
the induced strain and around the central part) is illus- 
trated by using Eq. (|2]). Figure [Sfa) shows a vector plot 
of the induced gauge fields where the length of the vec- 
tors and the colors denote the absolute value of A. The 
corresponding magnetic field is depicted in Fig.JHfb). No- 
tice that both the gauge field and the pseudo magnetic 
field exhibits an approximate six fold symmetry [l6|, [20| . 
Becasue of thermal fluctuations the symmetry is not ex- 
act. Notice that there is a little elongation along the sup- 
ported direction. We will discuss this symmetry in the 
next section. Notice that the induced magnetic fields are 
larger than those found for the step, trench and barrier. 
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FIG. 6: (Color online) The averaged (a) gauge field in and 
(c) the induced pseudo-magnetic field averaged over the y- 
direction for a graphene on top of a barrier as shown in Fig. [5] 
which has been supported from the longitudinal ends while it 
can freely move along the ^-direction. 



E. Gaussian bump/depression 

The Gaussian bump (protrusion) / depression [2l|, [22| is 
parameterized as 



hs(xi,yi) = ±/i exp(-p 2 2 /27 2 ), 



(8) 



where +ho(— ho)(=l nm) is the height (depth) of the 
Gaussian bump (depression) and pi 2 = x 



xf + yf is the 
radial distance of i th atom and 7 = 1 nm is the variance 
of the Gaussian. 

Since the optimum configuration of supported 
graphene over the Gaussian bump is similar to the one 
for a spherical bubble, we will not report them here. For 




FIG. 7: (Color online) The optimum configuration of arm- 
chair graphene, with two longitudinal ends supported in the 
x — y plane , on top of a bubble. The inset shows a different 
view indicating the elongation of the deformation of graphene 
along the x-direction, i.e. arm-chair direction. The colors 
indicate the size of strain. 
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FIG. 8: (Color online) (a) Vector plot of the gauge field and 
(b) the induced pseudo-magnetic field for arm-chair graphene 
placed over a spherical bubble. The obtained deformation is 
shown in Fig. [71 Graphene was supported from the longitu- 
dinal ends while it can vibrate along the z-direction. 



supported graphene over a Gaussian depression the op- 
timum configuration is not Gaussian (as was shown in 
Ref. [13). 

The results of the gauge fields and pseudo magnetic 
fields for a graphene membrane with Gaussian deforma- 
tion show a clear six fold symmetry [16, 20]. For a Gaus- 
sian deformation of the graphene membrane, Eq. flSJ, us- 
ing Eq. (|2j) and Eq. (0 we found 



A = 



2 7 4 



(cos(20),sin(20)/2), 



and for the corresponding pseudo-magnetic field 

h 2 (p,e) P 2 



B = V x A : 



■ sin(30) 



(9) 



(10) 



where x = pcos(#) and y = psin(#). The well known six 
fold symmetry is due to the dependence of B on sin(30). 
Our atomistic results for B (see Fig. |5fb)) are in good 
agreement with Eq. (fTD|) . 

It is surprising that we found a six fold symmetry for 
deformed graphene over a Gaussian bump but not for 
the Gaussian depression. This is due to the non-Gaussian 
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FIG. 9: (Color online) (a) Vector plot of the gauge field and 
(b) the induced pseudo-magnetic field for arm-chair graphene 
deposited over a Gaussian depression. Graphene has been 
supported from the longitudinal ends while it can vibrate 
along ^-direction. 



profile of GE on top of a Gaussian depression [15| . Break- 
ing the radial symmetry of graphene deformation reduces 
the symmetry in B (see Fig. [9](b)). This particular sym- 
metry affect also the energy eigne values and correspond- 
ing wave functions [l6| . 

In Fig. [9^a) the induced gauge fields for GE on top of 
a Gaussian depression is shown where the length of the 
vectors and the colors denote the absolute value of A. 
Both gauge fields and pseudo-magnetic fields are smaller 
than those found for the supported graphene over the 
bubble and the Gaussian bump. This is a consequence 
of the non-Gaussian profile for the optimum configura- 



tions which yields alteration in six fold symmetry in the 
induced pseudo-magnetic field. 



V. SUMMARY 

We investigated systematically the induced pseudo- 
magnetic field properties for graphene deposited on top 
of different nanostructured substrates by using molecular 
dynamics simulations at T=300K. The van der Waals 
interaction between the substrate and graphene was 
modeled by a Lennard- Jones potential. We found that 
the induced magnetic filed for graphene on top of a step 
consists of two magnetic barriers with different sign, 
while for a trench it forms two narrow magnetic barriers 
around the trench boundaries and one with opposite sign 
within the trench. The one directional substrate barrier, 
forms two sets of magnetic-barriers around the barrier 
wall. The magnetic filed for the Gaussian depression 
looses its six-fold symmetry (due to the non-Gaussian 
deformation of graphene) as compared to GE on top of 
the Gaussian bump/buble. The strain induced strong 
pseudo-magnetic fields. Controlling the pseudo-magnetic 
field is possible by controlling the substrate pattern and 
the size of the deformation. 
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